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SUMMARY 


The reciprocal equations relating the lift and moment acting on an 
airfoil due to harmonic oscillations in compressible flow to the indicial 
lift and moment functions for a pitching airfoil are determined. Two 
indicial functions are required to describe completely the lift or 
moment on the airfoil ; one function is used for describing the effects 
due to the angular position of the airfoil whereas the second function 
is used for describing the effects due to pitching velocity. 

Calculations for the indicial lift and moment ' functions are made 
for a pitching airfoil at a Mach number of 0.7 for which sufficient 
oscillatory coefficients were available . The growth of lift to the 
steady-state value is less rapid for compressible flow than for incom- 
pressible flow for both indicial functions of the pitching airfoil. 
Although the circulatory lift due to pitching velocity on an airfoil 
rotating about its three-quarter-chord point is nonexistent at a Mach 
number of 0, for a Mach number of 0.7 a time -dependent lift exists 
because of pitching velocity and approaches zero in approximately 
2 chords. 


INTRODUCTION 


In the study of transient flows, two types of airfoil motions have 
special significance - a harmonically oscillating airfoil and an air- 
foil experiencing a sudden change in angle of attack. The reciprocal 
relations between the indicial lift and the lift associated with a 
harmonically oscillating airfoil for both sinking and pitching motion 
are given in reference 1 for incompressible flow. In reference 2 the 
reciprocal relations were extended to apply to an airfoil suddenly 
acquiring a vertical velocity in a subsonic compressible flow, ahd the 




2 


NACA TN 2613 


indicial lift and moment functions were evaluated for a Mach number . 
of 0.7. The present paper is an extension of reference 2 in that the 
indicial lift and moment functions for a pitching airfoil are determined. 
The reciprocal equations for the case of a pitching airfoil are indi- 
cated, and the indicial lift and moment functions are computed for a 
Mach number of 0.7- 

In a recent paper by Lomax, Heaslet, and Sluder (reference 3)> a 
different method for determining the indicial lift and moment functions 
is given. Although parts of the indicial functions can be determined 
readily, the solution for the complete indicial functions is lengthy 
and tedious, and, consequently, numerical results are given only for a 
Mach number of 0.8 in reference 3- Part of the solution presented herein 
for a Mach number of 0-7j however, was determined by the method in refer- 
ence 3 and was compared with the solution determined in this paper. 


SYMBOLS 


s distance traveled, half- chords 

co angular frequency 

V • forward velocity of airfoil 

c chord 


x,z Cartesian coordinates 

k reduced- frequency parameter (coc/2V) 

L(s) lift per unit length of span, positive in downward 

v direction 


y»> 

M( s) 


M q (s) 


lift per unit length of span due to pitching velocity, 
positive in downward direction 

moment per unit length of span about quarter-chord 
point, positive when moment tends to depress 
trailing edge 

moment per unit length of span due to pitching velocity 
about the quarter-chord point, positive when moment 
tends to depress trailing edge 


P 


density 
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w 


norma l , component of perturbation velocity 


h 


amplitude of vertical displacement of airfoil, half- 
chords 


e 


angle of pitch, positive when trailing edge is lower 
than leading edge 


q 


ki(s) 


m L (s) 


kl q (s) 


rate of change of pitch of airfoil with respect to 
distance traveled in half-chords, positive when 
trailing edge is falling with reference to leading 
edge (dc/2V) 

indicial lift function for an airfoil experiencing a 
sudden change in vertical velocity 

indicial moment function for an airfoil experiencing 
a sudden change in vertical velocity where moment 
is taken about quarter-chord position 

indicial lift function for an airfoil experiencing a 
sudden change in pitching velocity about its 
leading edge 


mi (s) 


indicial moment function for an airfoil experiencing 
a sudden change in pitching velocity about its 
leading edge 


F c (k) + iG c (k) coefficient of complex compressible oscillatory lift 

derivative 

M(k) + iN(k) coefficient of complex compressible oscillatory 

moment derivative where moment is taken about 
quarter-chord point 

F c (k) + iG c (k) coefficient of complex compressible oscillatory lift 
^ ^ derivative due to pitching velocity only 


Mq(k) + N q (k) 


coefficient of complex compressible oscillatory moment 
derivative due to pitching velocity only 


f(k) = F c (k) - F c (») 


m(k) = M(k) - M(°°) 
fq(k) = F c (k) - F (») 
m q (k) = Mq(k) - Mq(») 
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Z 


l’ Z 2 


in-phase and out-of-phase lift coefficients, 

respectively, associated with translation of airfoil 






m 3 ,m 4 


in-phase and out -of -phase moment coefficients, 

respectively, about quarter-chord point associated 
with translation of airfoil 

in-phase and out -of -phase lift coefficients, 

respectively, associated with pitching motion of 
airfoil about its leading edge 

in-phase and out-of -phase moment coefficients, 

respectively, about quarter-chord point associated 
with pitching motion of airfoil about its leading 
edge 


METHOD OF ANALYSIS 


Consider the perturbation velocities on a two-dimensional wing 
(see the following sketch) for the case of combined vertical motion and 
pitching motion. For convenience the reference axis is taken as the 
leading edge and z and 6 are used to describe, respectively, the 
vertical and angular positions of the airfoil. 



The perturbation velocity (w measured positive upward, z positive 
downward) may be expressed as 


w(x,t) = -(TO + z + x0) 

If the airfoil is copsidered to have vertical motion only, then 
the perturbation velocity will be uniform across the chord and will be 
of intensity z. If the case of an airfoil pitching about its leading 
edge is considered, however, the perturbation velocity will be composed 
of two parts: a uniform part of intensity ■ TO and a linearly varying 
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part of intensity x0. When the lift on an airfoil having arbitrary 
motions is determined, therefore, these observations mean that only two 
types of perturbation velocities have to be considered: a uniform dis- 

tribution which applies either to vertical motion or to angular position 
alone and a linearly varying distribution which applies to angular 
velocity of the airfoil. 

When unit step motions of the airfoil are considered, it should be 
apparent that one indicial function is sufficient to define the lift 
on an airfoil suddenly acquiring a vertical velocity, since only the 
uniform perturbation velocity is involved. For the case of a pitching 
airfoil, however, two indicial functions are required: one corresponding 

to the uniform perturbation velocity associated with angular position 
and the other corresponding to the linear variation which is associated 
with angular velocity. 

The lift on an airfoil suddenly acquiring a vertical velocity is 
generally given in terms of an indicial lift function ki(s) by means 
of the equation 


L(s) = -JtpcV 2 H k^s) 


(la) 


where 


dh c_ * _ 2 

ds ~ 2V n ~ V 

This equation may be made to apply to an airfoil suddenly acquiring an 
angular position, because of the similarity that exists between the per- 
turbation velocities for the case of vertical motion and for the case of 
angular position alone. The lifts for these two cases therefore will 
be equal if the intensities of perturbation velocity are equal, that is, 
if z = V0. With this condition and equation (la), the lift on an air- 
foil following a sudden change in angular position is 


L(s) = -JtpcV 2 0kj_(s) (lb) 

For the pitching case then, the only indicial lift function that remains 
to be determined is the function associated with pitching' velocity of 
the airfoil. The equation for lift following a sudden change in angular 
velocity 0, where this angular velocity is, for convenience, taken 
about the le’ading edge, may be written similar to equation (la) or 
(lb) 

Lq( s) = -2JTPCV 2 qk]_^(s) 
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where q = ^0 and In (s) is the indlcial lift function associated 
a 2V q 

with a sudden acquisition of an angular velocity about the leading edge. 


The indicial lift and moment functions due to pitching velocity 
can he determined with the aid of the appropriate reciprocal equations 
and flutter coefficients in a manner similar to that shown for the 
sinking airfoil in reference 2. The available values for the lift and 
moment on the oscillating airfoil for pitching motion at subsonic Mach 
numbers, however, have been obtained numerically only for the total lift 
and moment; that is, the components due to angular position and pitching 
velocity have not been separated. Fortunately, the component due to 
angular position may be subtracted out of' the total oscillatory -coeffi- 
cients because of the similarity that exists between the lift on a 
sinking airfoil and the component of lift associated with angular posi- 
tion of a pitching airfoil. The available data for an airfoil oscil- 
lating harmonically in vertical motion therefore are used in the reduc- 
tion of the data for an airfoil oscillating in pitch to give the lift 
and moment associated with oscillatory pitching velocity alone. These 
data can then be used to determine the indicial functions associated 
with pitching velocity. 


Reciprocal Equations 

The reciprocal relations given by equations (8a) and (8b) in refer- 
ence 2 between the compressible lift on an airfoil due to harmonic 
oscillations and the indicial lift on an airfoil experiencing a change 
in angle of attack due to a sudden acquisition of vertical velocity are 
as follows: 


and 


ki(s) 



F c (k)sin ks 
_ 


dk 


(s > 0) 


(2a) 


k^s) 


F c (0) 



G c (k)cos ks 
k 


dk 


(s > 0) 


(2b) 


where F c (k) and G c (k) are related, respectively, to the in-phase and 
out-of -phase lift components on an oscillating airfoil and are defined 
by the following equation: 


L(s) = -itpcV 2 e iks (-ikh) 


r. 


l f c (s) 


iG r (k) 


) 


(3) 
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For the case of a pitching airfoil the component of lift associated 
with angular position and its corresponding indicial lift function are 
given hy equations (2) and (3)* where the term 


-ikh = 6 




a relationship which follows from the fact that the perturbation veloci- 
ties for vertical motion are equal to the velocities due to angular 
position alone when z = VS. 


In the interpretation of the expression for lift given by equa- 
tion (3), it should be noted that the expression for the motion of the 
airfoil is in complex form. If, for example, the actual angular posi- 
tion of the airfoil is denoted by 


I.P.Se 


iks 


then the lift is 


I.P.-itpcV 2 0e iks F c (k) + iG c (k) 
_ — 


or 


-npc 


V 2 © 


F c (kr + G c (kr 


01V2 


J 


sin 


ks +■ tan 


- 1 . 


G c (k) 

FjkJ 


Thus, the lift and the angular displacement have the same frequency and 
both their magnitude and phase are functions of k. 

The reciprocal equations for the component of lift due to pitching 
velocity may be expressed in terms of a harmonically oscillating airfoil 
and the indicial lift function due to a sudden change in pitching 
velocity kq (s) as follows: 


P 00 F (k) sin ks 




(s >0) (5a) 


and 


kq (s) = F c (0) + | / 

q q * J 0 


2 / c q 


G c -(k) cos ks 


dk (s > 0) (5b) 
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where F C g(k) and G C g(k) are related, respectively, to the in-phase 

and out-of -phase lift components on an oscillating airfoil due to 
pitching velocity about the leading edge and are defined by the equation 


L q (s) = -irpcV 2 e ik8 (2ik0)[F c ^(k) + iG Cg (k)j (6) 

In equations (5a) and (5b) the functions F c ^(k) and G C g(k)/k must 
be continuous and finite in the interval from k = 0 to k = <». 

The available data for evaluating the F c (k), G c (k), F c (k), 

H 

and G c ^(k) functions have been given in various forms, beginning with 

the work of Possio. One form in which all of the available data can be 
put is as follows: 

Translatory motion: 


L(s) = JtpcV 2 e ikB + iZ 2 ) (7a) 

■ M(s) = Jtpc^V 2 e 1 ^ 8 + iM^ (7h) 

Pitching motion (including contribution due to pitching velocity): 

L(s) = -npcV 2 e iks e(z 3 + iZ^) (8a) 

» 

M(s) = -rtpc^e 11518 0 (m 3 + iM]^ (8b) 


where the moment is taken about the quarter-chord point and the axis of 
rotation for the harmonically pitching airfoil is located at the leading 
edge. 


In order to make use of the reciprocal equations (5a) and (5b) to 
determine the indicial lift function due to pitching velocity kig(s), 

the component of the flutter coefficients due to pitching velocity must 
first be separated from the flutter coefficients given by equations (8a) 
and (8b) and then converted to the form given by equation (6). If the 
expression for h determined by equation (4) is substituted into 
equation (7a) and if the resulting expression is subtracted from equa- 
tion (8a), the following expression for the lift can be obtained for the 
component due to pitching velocity on a harmonically oscillating airfoil: 
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V s ) 


-itpcV 2 e iks 




(9) 


Comparison of equations (6) and (9) lead? to the following expressions 


for the functions 
coefficients Z 


F c .(k) and 




Z3, and 


Gc q (fc) 


in terms of the oscillatory 




and 


if ^1 
F c = Skl z i<- + Sk, 

n ' 


(10a) 


• 1 / 2ip' 

“-a( z 3 -m 


(lot) 


Substitution of equations (10) into equations (5) then allows for the 
determination of the indicial function kq (s) for pitching velocity. 


The procedure shown for determining’ the indicial lift function due 
to pitching velocity from the oscillatory coefficients may be applied 
to determine the indicial moment function due to pitching velocity mq (s) 

defined as follows: 


M(s) = 2rtpc 2 V 2 <5a 1 (s) (11) 

' . q 

where q is the magnitude of the pitching velocity. The reciprocal 
equations can be written for the moment in terms of the oscillatory 
coefficients and the indicial moment function mi (s) as follows: 


2 r ro M q (k) sin ks 

mq^s) = nj — 5 ( s > °) ( 12a ) 


and 


m l q ( s ) = Mq(0) + | J 


2 r°° Kq(k) cos ks 


dk (s >0) (12b) 


where M^(k) and Nq(k) are, respectively, the in-phase and out-of- 
phase moment coefficients on an oscillating airfoil due to the angular 
pitching velocity about the leading edge and are defined by the equation 
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Mq(s) = Jtpc 2 V 2 e iks (2ik0) jM q (k) + iN q (k)j (13) 

In equations (12a) and (12b) the functions M^(k) and Nq(k)/k must 
be continuous and finite in the interval from k = 0 to k = co . 


In a manner similar to tbat shown for the lift case, expressions 
for M n (k) and N a (k) may be obtained in terms of the oscillatory 

SI 

coefficients M-^, Mg, M^, and M^ as follows: 

V k) - -§s( m 4 + 1) (l4a) 

and 


V k > 




(ltb) 


Summary of Available Flutter Coefficients and Related 
Data for a Mach Number of 0.7 

The two-dimensional compressible flutter coefficients for the real 
and imaginary parts of the lift and moment for sinking and pitching 
motion are given in table I for a Mach number of 0-7- The results given 
in table I were taken from three sources (references 4 to 6) and have 
been converted to the form given by equations (7) and (8) for the lift 
and moment. The range of reduced frequencies taken from each source is 
indicated in table I. 

Expressions for the end points F c (0) and F c (°°) can be deter- 
mined independently of the flutter coefficients. The value of F c (0) 
which represents the steady-state lift is determined by the Prandtl- 

Glauert factor — — (where M refers to the Mach number of the 

free stream) . The value of F c (<») which was shown in reference 2 to 
correspond to the value of k^( s) at s = 0 can be determined by the 
following equation given in reference 2 for any Mach number: 
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The corresponding expression for the moment M(co) is noted in refer- 
ence 2 as 


M(°°) 


1 

2M 


(16) 


For the case of a pitching airfoil the expression for the end points 
for the F c (k) and Mq(k) functions may he obtained in a manner similar 

to that shown for the case of sinking airfoil. For an airfoil exper- 
iencing a sudden change in pitching velocity about its leading edge, 
the steady-state value of the kp^(s) function (which can be shown to 

correspond to F c (0)) can be determined from incompressible flow 

Si 

together with a correction given by the Prandtl-Glauert factor. From the 
work of either Theodorsen or Wagner (reference 8 or 9) the incompres- 
sible steady-state value can be determined, and, together with the 
Prandtl-Glauert factor, the resulting expression for F C( ^(o) for 

subsonic flow is as follows: 


F c (°) = kj. (») = 1 






Similarly, for the moment the expression for M q (o) is 


/X , . 11 

M q(°) = % (~) = -35 — 

q l ±D \J1 - M 2 


(17) 


(18) 



values of the k-i (s) and m n (s) functions, respectively. The same 

1 J -q. 


method of analysis for obtaining the starting values of the kp(s) 
and m^(s) functions given in reference 9 for a sinking airfoil can be 
applied to the case of a pitching airfoil to obtain the starting values 
for the kn (s) and mn (s) functions if a substitution of the boundary 

q q 

conditions due to pitching velocity is made in lieu of sinking velocity 
as indicated by equation (l). If the perturbation velocity component 
due to a sudden change in pitching velocity is substituted for the 
perturbation velocity due to sinking motion in the analysis given in 
reference 9 f or determining kj_(s) and mq(s) at s = 0, the following 
expressions are obtained for the starting values of the kp^(s) and 

mi^(s) functions 

kp (0) = 

x q JtM 


( 19 ) 
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and 




( 20 ) 


The values of 


kp^(o) and mp^O), however, correspond, respectively. 


M q (k) 


at 


to the values of the flutter coefficients F c ^(k) and 

infinite frequency. The proof for this correspondence is identical to 
the one shown for the sinking airfoil given in appendix A of refer- 
ence 2. Thus equations (l9) and (20) denote, respectively, the values 


of F„ (°°) and 


M q (~) 


and 



( 21 ) 


M q (°°) =- 


5_ 1_ 

12 jtM 


( 22 ) 


RESULTS AND DISCUSSION 

Numerical Solution of the Reciprocal Equations 
for a Mach Number of 0-7 


As noted in reference 2, the indicial lift and moment functions 
determined by the reciprocal equations containing the in-phase lift or 
moment functions (F c (k) or M(k)) provided a more reliable and simpler 
solution than the functions determined by the reciprocal equations con- 
taining the out-of -phase components (Gc(k) and N(k)). Consequently, 
the indicial lift and moment functions for the pitching-velocity case 
were determined numerically by using equations (5a) and (12a) in a man- 
ner similar to that shown for the case of the sinking airfoil. 

In figure 1 a plot is shown of the complex oscillatory lift func- 
tions due to pitching velocity at a Mach number of 0.7 which were 
evaluated from the flutter coefficients by use of equations (10a) and 
( 10b) . If 

-?= q < k ’ -%<"> 


(23) 
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is substituted, into equation (5a) , the following form which can be 
evaluated graphically more readily is obtained: 




'“f n (k) sin ks 


dk 


(24) 


A plot of the function fq(k) is shown in figure 2. The integrand in 
equation (24) was graphically evaluated for several values of the parame- 
ter s. .In figure 3 , plots of the indicial lift function kq^s) are 

shown for Mach numbers of 0.7 and 0. 'The indicial lift function kq(s) 
due to angular position alone (corresponding to the case of a sinking 
airfoil) is plotted in figure 4° as determined from reference 2 for a 
Mach number of 0.7. Also shown in this figure are the beginning posi- 
tion of the curve for a Mach number of 0.7 as determined by the equa- 
tions given in reference 3 , the curve for a Mach number of 0.8 given 
in reference 3, and the curve for a Mach number of 0. Comparison of 
the kq(s) , and kq^(s) functions at Mach numbers of 0 and 0.7 indi- 
cates that the growth of lift to the steady state for a pitching airfoil 
is less rapid in subsonic compressible flow than for incompressible flow. 

The indicial moment function mig(s) due to a sudden change in 

pitching velocity may be obtained in a manner similar to that shown 
for determining the indicial lift function kq q (s). In figure 5, plots 

are shown of the complex moment function due to pitching velocity for a 
Mach number of 0.7 determined from the flutter coefficients Mq, M2, 

Mg, and M^ by equations (l4) . The asymptotic value for M q (k) is 
determined by equation (22) . Because of the nature of the curves at 
the higher reduced frequencies, the curves are represented by dashed 
lines to indicate a degree of unreliability. Reference 2, however, 
indicates that errors in the curves at the higher reduced frequencies 
would have little effect on the indicial function. 

In a manner similar to that shown for the lift case, the M q (k) 
function was transformed by the following equation 


m q (k) = M q (k) - M q («) (25) 


in order to facilitate the numerical determination of the indicial 
moment functiop mq (s). A plot of m_(k) is shown in figure 6. The 


function 


mi q (s) 


was evaluated graphically by substitution of equa- 


tion (25) into equation (12a) for several values of the parameters and 
is plotted in figure 7 together with the solution for a Mach number 
of 0. In figure 8 a plot of the 'indicial moment function due to 
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angle-of-attack change (no pitching velocity) mp(s) is shown for Mach 
numbers of 0.7 and 0.8 determined, respectively, from references 2 
and 3 together with the solution for a Mach number of 0. 

References 8 and 9 show that the indicial lift due to pitching 
velocity on an airfoil pitching about its three -quarter- chord position 
is impulsive at s = 0 and is zero thereafter for incompressible fieri'?. 
For compressible flow the indicial lift, however, is initially finite 
and is time-dependent thereafter, as is shown subsequently. When the 
results presented in figures 3> 4, 7 , and 8 are- used, the indicial lift 
and moment functions due to pitching velocity for an airfoil pitching ■ 
about any axis can be obtained readily by the following general trans- 
formations : 


( kl q) x ( s ) = k l q ( s ) ' f k l( s ) 

and 

(%)x (s) = " c mi(s) 


(26a) 


(26b) 


where x denotes the location of the axis of rotation measured posi- 
tively to the right from the leading edge of the airfoil (see sketch in 
section entitled "Method of Analysis") . For the case of an airfoil 
pitching About its three-quarter-chord location, equations (26a) and 
(26b) become 



= k-L (s) 

q 




(s) 


and 


(27a) 



m l q ( s) 



(27b) 


The functions 


N 


3c 

4 


( s) and 


Nso 


(s) were determined for a 


Mach number of 0-7 from equations (27a) and (27b) and from the results 
plotted in figures 3> 4, J. and 8. Plots of the functions /in ) (s) 

V v 3c 

4 
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and 


( miq )i£ 


(s) are shewn in figures 9 and 10, respectively, together 


with the indicial functions for the case for a Mach number of 0. 

Also included in these figures are the solutions for a Mach number 
of 0.8 and part of -the solution for a Mach number of 0.7, both deter- 
mined from reference 3* A comparison of the ^kp j (s) function for 


subsonic compressible and incompressible flow indicates that, although 
this function is impulsive at s = 0 and zero for s > 0 for incom- 
pressible flow, the function is finite at s = 0 and decays to zero 
very rapidly in the compressible case. If for compressible flow the 
circulation is assumed to be zero for s > 0 (as is the case for this 
component at M = 0) , the time-dependent function which is present in 
compressible flow may be attributed to the time -dependent apparent-mass 
effects. Also, the part of the kp (s) function shown in figure 4 

for s greater than approximately 4 may be associated only with the 
lift due to circulation. ■ Comparison of this part of the curve with the 
curve for a Mach number of 0 in the region s > b indicates that the 
growth of lift to the steady state is relatively less rapid for com- 
pressible flow than for incompressible flow. 


Approximation of Indicial Lift and Moment Functions 
by Analytic Expressions 

Since the exponential function has a simple operational equivalent 
and it has been found convenient to approximate the kp(s) function at 
a Mach number of 0 (see reference 10) , a limited series of such func- 
tions were chosen to approximate the indicial lift and moment functions 
for sinking and pitching motion at a Mach number of 0.7. The functions 
were found to fit these curves quite well and are: 

k^s) = 1.4(1 - 0.361k- 0 - 0536s - 0.4O5e-°- 35 ' 7B + 0.4l9e-°- 902s ) ( 28 ) 

fe) <») - -0.083e-°- 8006 - O^e" 1 ' 51556 + 0.l49e' 2 - Ws (29) 

T 

mp(s) = -0.2425e" 0 ’ 9 ^ 8 + 0.084e -0 ’ 668s - 0.0 69e -0 ‘^ 38s (30) 
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-0.0875 (] 


1 + 0.1l4le 


-0.1865s 


1.233e 


-1.141s 


+ 0.3337e 


- 4 .o 4 s\ 


( 31 ) 


These analytic approximations given "by equations (28), (29), (30), 
and (31) are plotted in figures 4, 9, 8, and 10, respectively. Com- 
parison of the approximation with the actual curves in these figures 
indicates the order of good agreement reached. 

The corresponding approximate expressions for the harmonically 
oscillating airfoil can he found from the reciprocal form of equa- 
tions (2) and (5) where the harmonically oscillating functions are 
expressed in terms of their corresponding indicial function. For the 
case of the sinking airfoi-1, the reciprocal form of equation (2) is 


F c (k) + iG c (k) = fk 



ds 


(32) 


If ik is considered as the operator in the Laplace transformation, . 
then equation (32) is simply ik times the Laplace transformation 
of kp(s). Therefore F c (k) and G c (k) are, respectively, the real 
and imaginary parts of ik times the Laplace transform of kj( s) . 
Substitution of equation (28) into equation (32) leads to the following 
expressions for F c (k) and G c (k) : 


F_(k) = 1.4 


1 - 


■364k 2 


0 . 405 ^ 


. 419k 2 


(0:0536) 2 + k 2 (0.357) 2 + k 2 (0.902) 2 + k 2 


> ( 33 ) 


G c (k) 


1.4 


(0.364) (o.Q536)k 
(0.0536) 2 + k 2 


(0.405) (0.357)k 

(0.357) 2 + k 2 


(Q.4i9)(o.9Q2)k 

(0.902) 2 + k 2 J 

J 


In a manner similar to that shown for the lift case, the following 
approximate expressions for M(k) and W(k) can he obtained: 
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IT 


M (k) = -0.2425k 2 _ + 0.084k 2 


0.069k 2 




(0.974) 2 + k 2 ( 0 . 668 )* + k 2 (0.438) 2 + k' 


2 , ,2 


,2 . ,2 


N ( k ) = ( -0.2425) (0. 974) k + (0.084) (0.668)k _ (0.069) (o.438)k 
(0.974) 2 + k 2 (0.668) 2 + k 2 (0.438) 2 + k 2 


(34) 


j 


The expressions for the k ) and ^^^(k) lift functions 

T T 

due to pitching velocity are 


fa \ ( k ) = -0.083k 2 0.293k 2 + 0.149k 2 

V q /3c (0.800) 2 + k 2 (I.565) 2 + k 2 (2.44) 2 + k 2 
4 


/ G \ ( k ) = (-0.083) (Q -8oo) k _ (0.293) ( 1^569) k + (0.149) (2. 44)k 
\ ° q '3c (0.800) 2 + k 2 (1.565) 2 + k 2 (2.44) 2 + k 2 


> (35) 


and the expressions for the (k) and ^Nqj^k) 

T 4 

are 


moment functions 


^M q j 3c (k) = -0.0875 
4 


1 + 


o.n4ik 2 


1.233k 


(0.1865) 2 + k 2 (l.l4l) 2 + k 2 (4.04) 2 + k : 


0.3337k 


2 . ,2 


v2 . ,2 


W 3c (k) 

T* 


= -0.0875 


(0.1141) (0. 1869) k (1.233) (l.l4l)k 


j_ (0.1865) 2 + k 2 (l.l4i) 2 + k 2 


(0.3337) (4. o4)k 
(4.o4) 2 + k 2 


(36) 
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CONCLUDING REMARKS 


The indicial lift and moment functions due to a sudden change in 
pitching velocity have been obtained from the available data on a 
harmonically pitching airfoil at a Mach number of 0.7 by the use of 
reciprocal relations. Comparison of the results obtained for the case 
of a Mach number of 0.7 with incompressible flow indicates that the 
growth of lift to the steady' state appears to be less rapid for com- 
pressible flow than for incompressible flow for a pitching airfoil. 
Although the circulatory lift component due to pitching velocity for 
an airfoil rotating about its three -quarter- chord point was zero for 
incompressible flow, a time-dependent component of lift was found to 
exist for compressible flow and decayed to a negligible value in 
approximately 2 chords. 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va., October 30 , 1953 
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TABLE I 


ro 

o 


LIFT AND MOMENT OSCILLATORY COEFFICIENTS FOR SINKING AND 
PITCHING MOTION AT A MACH NUMBER OF 0.7 


Reference 

k 

mm 

z 2 

z 3 

Z4 

M 1 

Mg 

m 3 

mm 


0 

0 

0 

1.4003 

0 

0 

0 

0 

0 


.02 

.00711 

.05177 

1.2985 

-.13891 

-.00026 

.00003 

.000574 

.010002 


.04 

.02006 

.09524 

1.2059 

- . 17914 

-.00097 

.00016 

.001452 

.019332 

- 

.06 

.03399 

.1324 

1.1304 

-.18304 

-.00208 

.00042 

.002319 

.028182 


.08 

.04722 

.1651 

1.0698 

-.17128 

-.00355 

.00083 

.003090 

.036758 

4 

.10 

.0589 

.1941 

1.0216 

-.14992 

-.0053 

.00126 

.003675 

.045114 


.20 

,0945 

.3186 

.8869 

-.003247 

-.0186 

.00574 

.004301 

..086501 


•30 

.0991 

.4332 

.8364 

.14614 

-.0393 

,01346 

.002632 

•12953 


.4o 

.0827 

• 5523 

.8199 

.28495 

-.0672 

.02596 

.000512 

• 17541 


.50 

.0538 

.6820 

.8239 

.4135 

-.1014 

.04482 

.000436 

.2242 


.60 

.0187 

.8229 

• 8399 

.5294 

-.1402 

.07175 

.000573 

.27498 


.70 

-.0154 

• 9751 

.8669 

.6355 

-.1813 

.1090 

.010257 

.32692 

5 

.80 

-.0449 

1.1316 

.8934 

.7191 

-.2178 

.1578 

.021912 

•37005 


1.00 

-.0874 

1.4482 

• 9391 

.8607 

-‘.2798 

.2665 

.059166 

.45875 


1.50 

-.39295 

2.311 

.8625 

1.3512 

-.44oo 

.4896 

-.001739 

.58628 

6 

2.00 

-.6521 

3.306 

.8110 

1.8258 

-.5740 

•77185 

-.062784 

.73568 


2.50 

-.6083 

4.241 

.8566 

2.0913 

-.4547 

1.218 

.090049 

•93504 




NACA TN 2613 



































Figure 3 .- Comparison of indicial lift functions due to a sudden change 

in pitching velocity for an airfoil rotating about its leading ddge ro 

at M 0 and M = 0.7. w 
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Figure 4.- Comparison of indicial lift functions due to a sudden change 
in angle of attack (without pitching motion) at M = 0, M = 0.7, and 
M = 0.8. 








Figure 5.- Oscillatory moment function due to pitching velocity for an 
airfoil pitching about its leading edge at M = 0.7, (Moment taken 
about quarter-chord point. ) 









Asymptotic value 
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M=0- 


s, hair-chords 


Figure 7*~ Comparison of indicial moment functions due to a sudden change 
in pitching velocity for an airfoil rotating about its leading edge 
at M >= 0 and M = 0.7* (Moment taken about quarter-chord poinlT. ) 
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* Figure 9 .- Comparison of indie ial lift functions due to a sudden change 
in pitching velocity for an airfoil rotating about its three-quarter- 
chord point at M = 0, M = 0.7, and M = 0.8. 
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! Figure 10.- Comparison of indicial moment functions due to a sudden change 

j g in pitching velocity for an airfoil rotating about its three-quarter- 

I > chord point of M » 0, M = 0.7, and M = 0.8. (Moment tafeen about 

J | quarter-chord point. ) 

! f 

I I 


i 


MACA TN 2613 




